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Abstract. - The around-the-mean-field version of the LDA+U method is applied to investigate 
electron correlation effects in 5-Pu. It yields a non-magnetic ground state of 8— Pu, and provides 
a good agreement with experimental equilibrium volume, bulk modulus and explains important 
features of the photoelectron spectra. 



Plutonium is probably the most intricate element from the point of view of condensed 
matter physics. It exhibits six allotropic modifications at ambient pressure, some of them 
of very low symmetry (monoclinic) , and there is little doubt that the anomalous behaviour 
is related to the 5/ electronic states, being at the borderline between the localized, non- 
bonding, behaviour and the bonding situation of electronic bands. One can preconceive that 
at the cross-over regime, several states with very different degree of 5/ derealization can be 
nearly degenerated in energy. 

The most thoroughly studied phases are a— Pu (monoclinic) and 6— Pu (fee). The latter 
is stable between 592 and 724 K, but can be stabilized down to T = K by various dopants. 
This phase has the largest volume (by 20% higher than a— Pu) (for overview see Ref. [1] and 
references therein). 

The atomic volume is an important indicator of the situation of the 5/-electronic states. 
Withdrawing of the 5/ states from the bonding and confining them in the ionic core leads to a 
significant volume expansion. On the plot of atomic volumes of elements, Pu (a) represents a 
continuation of light actinides, with the decreasing branch resembling the parabolic behavior 
of transition metals. On the other hand, heavy actinides, starting with Am (Z = 95), display 
higher volume, following a weakly decreasing volume of lanthanides, characterized by non- 
bonding 4/ states. 6— Pu, being half way between the volume of a— Pu and Am, represents 
therefore generally the cross-over regime, where electron-electron correlations play a prominent 
role. 

Ab-initio electron energy calculations based on the Density Functional Theory (DFT) in 
the Local Density (LDA) or Generalized Gradient (GGA) approximations account generally 
well for basic properties of metallic systems. Numerous variants of this successful paradigma 
were applied to Pu phases. The most conspicuous failure is the case of 8— Pu calculations which 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



inevitably lead to magnetic ordering. The fact that the lattice expansion due to magnetism 
yields approximately correct value of the volume, and intrasingence of magnetic order within 
the DFT theory both within the LDA and GGA approximations, led to speculations about 
magnetic ordering of 5— Pu. But it contradicts experimental findings (magnetic susceptibility 
has a character of weak Pauli paramagnet [2], paramagnetic state is also evidenced by 27 Al 
NMR [3], 69 Ga NMR [4] and neutron scattering [5]). Since its first occurrence [6], new claims 
of magnetic Pu keep filling research journals; there is already a debate about the detailed 
type of order, and recently even other Pu phases were assumed as magnetic [7]. One common 
feature of light actinides is that orbital moments, appearing due to the strong spin-orbit 
interaction in the case of spin polarization, are oriented antiparallel to spin moments. LDA 
or GGA calculations of S— Pu (with approximately five 5/ electrons) lead to the spin moment 
larger than orbital moment (which is opposite to what the Hund's rules assume for the 5/ 5 
ionic state), but the existence of "exact" accidental complete cancellation can be doubted. 
The systematic cancellation appears only for the 5/ 6 configuration of Am (both for the LS 
and j j coupling) . 

Pu magnetism as an undesirable artefact tends to appear also in other approaches, like 
the mixed-scheme model, based on minimization of energy for the localized 5/ 4 configuration, 
while the remaining 5/ electrons form a band [8], while the Dynamical Mean Field Theory 
(DMFT) deduces that magnetic order is washed out by dynamical fluctuations of moments [9] . 
As yet, the challenging question and the real problem of Pu remains, which is why is there no 
magnetism in the <5-Pu phase? 

Here, we address this problem making use of a well known L(S)DA+U method. Generally, 
the LDA+U calculations account for the on-site correlations between the / electrons in a more 
realistic way than the LSDA. First applied to 5— Pu by Bouchet et al. [10], they also lead to 
a magnetic solution, similar to that reported by Savrasov and Kotliar [11]. In a contrary 
to [10, 11], we apply to S— Pu a different version of the LDA+U method which is based on 
the original LDA+U total energy functional of Ref. [12]. We show that when the LDA+U of 
Rcf. [12] is reformulated in a spin and orbital rotationally invariant form, it yields a basically 
non-magnetic S— Pu with S — > and L — > 0. 

The correlated band theory L(S)DA+U method consists of the local spin-density approx- 
imation (LSDA) augmented by a correcting energy of a multiband Hubbard type E ee and a 
"double-counting" subtraction term Ed c which accounts approximately for an electron-electron 
interaction energy already included in LSDA. It is well known that the form of "Coulomb-?/" 
correction to the LDA is not uniquely defined [13]. The most commonly used is the version of 
LDA+U total energy functional in a so-called "fully localized" limit (FLL) [11,14], in which 
the double-counting term Ed c is taken to satisfy an atomic-like limit of the LDA total energy. 

Another and historically the first LDA+U functional is often called as an "around-mean- 
field" (AMF) limit of the LDA+U. In this AMF-LDA+U limit [12, 15] the interaction energy 
takes the form 

E^ F = \ E <Ki,7 2 [<7i,73|V ee | 7 2,74>- (7i, 7 3 |V ee | 7 4, 72)] <5% 3 , 74 , (1) 

71,72,73,74 

where V ee is an effective on-site "Coulomb-?/" interaction, the combined spin-orbital index 
7 = (ma) is composed from the angular m and spin index a, and 

1 1 

^71,72 — ^71,72 — n 1< -*7l,72 1 n 2/ + 1 ^ ' n ma,ma i (2) 
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where n. 7li72 is the on-site /-occupation matrix in the spin-orbital space which has to be 
defined with respect to the chosen localized orbital basis set, and n a is an average spin-orbital 
occupation. The double-counting correction is set equal to zero, £^, MF = 0. 

We note that the essential feature of the Eq. is the presence of spin-off-diagonal 
elements of the on-site occupation matrix n 7l72 = n miai ^ m2a2 which are in general non-zero 
in the presence of the spin-orbit coupling (SOC). 

As it was shown recently in Ref. [13] , the AMF-LDA+U formulation is most appropriate for 
metallic phases in which the Coulomb-C/ is comparable to the bandwith W . On the other hand, 
the FLL-LDA+U model is more suitable for insulators, with Coulomb-C/ much larger than W . 
Here we have extended the AMF-LDA+U for the general spin and orbital rotationally invariant 
case including the SOC, and implemented it in the full potential linearized augmented plane 
wave (FP-LAPW) basis in a way similar to a previous relativistic FLL-LDA+U FP-LAPW 
implementation [16]. 

We calculated the electronic and magnetic structure of 5-Pu using three different ap- 
proaches: (i) LSDA, (ii) FLL-LDA+U, and (iii) AMF-LDA+U. In all these calculations, 
we assumed a tetragonal unit cell with two Pu atoms in order to accommodate the anti- 
ferromagnetic (AFM) order, as it is well known that both the LSDA [17] and GGA [18] yield 
the AFM as the ground state of Pu. 

Table I - Ground state spin Ms, orbital Ml, and total Mj = Ms + Ml magnetic moments (in ^b) 
in 5-Pu calculated for antiferromagnetic configuration at the experimental lattice parameter a — 8. 760 
a.u. using the LSDA, FLL-, and AMF-LSDA + U (U = 4 eV) models. Also given are the equilibrium 
volume V eq (in (a.u.) 3 ) and bulk modulus B. We note that at the equilibrium the AFM-LSDA is lower 
in the total energy than non-magnetic LDA by 12.2 mRy/atom. 



Model 


M s 


M L 


Mj 


v eq 


B (kbar) 


LSDA 


4.357 


-2.020 


2.337 


136.8 


761 


FLL LSDA+U 


3.272 


-3.802 


-0.530 


187.9 


675 


AMF LSDA+U 


~ 


~ 


~ 


181.5 


314 


Experiment 


N/A 


N/A 





168 


299 



The spin M$, orbital Ml and total Mj magnetic moments on Pu atom calculated within 
the LSDA (see Table I) are in a very good agreement with previous LSDA calculations [17]. 
Starting from the LSDA calculated charge and spin densities and on-site spin and orbital 
occupations, we performed the LDA+U calculations. We choose U — 4 eV and exchange 
interaction J = 0.7 eV in the range of commonly accepted values for Pu [11]. The FLL- 
LDA+U model yields a substantial (approx. 1 /j,b) reduction of Ms and an increase of \Mt\ 
by approximately 2 fig, resulting in a substantial reduction of the (see Table I), while 

keeping non-zero local magnetic moment on Pu atom [19]. The Coulomb-C/ has a strong effect 
on both the spin and orbital polarization, but it changes very little the total /-state occupation 
nt - from approximately 5.06 in LSDA to 4.95 in FLL-LDA+U. With the increase of U from 
4 eV to 5 eV, the Ms increased to 3.42 fi B and Ml decreased to —3.97 lib, yielding a slight 
increase of |Mj| to 0.55 fig and a decrease of nf to 4.93. Thus, in spite of major differences in 
the values of the Ms, Ml, and Mj, both the LSDA and FLL-LDA+U calculations describe 
5-Pu as a 5/ 5 state. 

Next, we turn to the salient aspect of our investigation, the AMF-LDA+U (U = 4 eV, J = 
0.7 eV) calculations. Starting from strongly spin-polarized LSDA charge and spin densities 
and 5/ manifold spin and orbital occupations and without any constraint, the calculations 
converged to the almost zero magnetic moment with remaining Ms and \M~l\ less than 0.01 
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Fig. 1 - Spin Ms, orbital Ml, and total Mj magnetic moments dependence on the Coulomb-(7 as 
calculated in AMF-LDA+U model for the experimental lattice parameter a = 8.760 a.u. 

Fig. 2 - Total energy as a function of the volume per Pu atom calculated within the AMF-LSDA+U 
((7 = 4, J — 0.7 eV) and LSDA in the antiferromagnetic model. 

Hb- We also performed the calculations starting from a different FLL-LDA+U ground state 
and obtained essentially the same results. We conclude that the AMF-LDA+U yields fun- 
damentally non-magnetic <5-Pu, in accordance with experimental observations. Importantly, 
the 5/ occupation nf is increased substantially from rif w 5 in LSDA and FLL-LDA+U to 
fif = 5.44 (see Table I) meaning that there is a substantial deviation from the 5/ 5 ionic state. 

The dependence of the spin Ms, orbital Ml, and total Mj calculated within the AMF- 
LDA+U on the value of Coulomb-C/ is shown in Fig. 1. For small values of U (~ J), <5-Pu is 
magnetic with sizeable Ms and Ml moments that almost cancel each other. As the U value 
is increased to 1.5 eV, the local moments Ms and Ml disappear and 5-Pu is non-magnetic 
for realistic values of the Coulomb-l/ (from 3 to 5 eV). It is interesting to note that rif also 
depends on U and it increases from n/=5.1 (17=0.5 eV) to n/=5.5 ([7=5 eV). 

The total energy dependence on the volume per Pu atom calculated within the AMF- 
LDA+U (f/ = 4 eV) is shown in Fig. 2 and compared with the results of the AFM LSDA 
calculations. First, we mention that the AFM LSDA yields the equilibrium volume V eq 136.8 
(a.u.) 3 per Pu and bulk modulus B = 761 kbar (see Table I), in a good agreement with 
the LMTO results of Ref. [17]. The LSDA V eq is close to a-Pu value 134 (a.u.) 3 which 
is approximately by 20 % smaller than the experimental volume 168 (a.u.) 3 of <5-Pu, while 
the bulk modulus is almost twice bigger than its experimental value 299 kbar [20] . As it was 
already shown in [10,11], the FLL-LDA+U model gives the V eq close to the experimental <5-Pu 
volume and it slightly lowers the bulk modulus (see Table I). Finally, the AMF-LDA+U yields 
V eq = 181.5 eV and B — 314 kbar (see Table I and Fig. 1) which are in a very good agreement 
with experimental data. As yet, the Coulomb-[7 is treated as an adjustable parameter and 
V eq can be fitted further to the experiment by a slight decrease of U. However, it is clear that 
the conventional value [7 = 4 eV works already quite well. Thus, the AMF-LDA+U is shown 
to give the best agreement with experiment both for magnetic and structural properties of 
the (5-Pu. 

In Fig. 3(a) we show the total (TDOS, per unit cell) and /-states projected densities of 
states (fDOS) resulting from the <5-Pu bulk AMF-LDA+U calculations. In order to clarify the 
character of the 5/ manifold resulting from the AMF-LDA+U calculations, it is convenient 
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Fig. 3 - a) Total (TDOS) and f-projected (fDOS) DOS as the results of AMF-LSDA+U=4 eV calcula- 
tions for the bulk <5-Pu; b) <5— Pu j z — resolved /DOS from AMF-LDA+U calculations; c) experimental 
UPS spectra for Pu multilayers and the bulk taken from Ref. [21,22]; d) 1ML, 3ML, and 5ML total 
DOS calculated for free-standing Pu multilayers making use of AMF-LDA+U model. 



to look at the j z = m; + m s /DOS shown in Fig. 3(b), as the z— projection of the total 
moment j z is the only remaining quantum number in the presence of the SOC and magnetic 
polarization. It shows six filled j — 5/2 /-states and the j = 7/2 /-states that are split 
away by approximately 4 — 5 eV. As the LDA+U model is not based on any kind of atomic 
coupling scheme (LS or jj) it rather determines the set of single-particle orbitals that minimize 
variationally the total energy. The LDA+U result can be interpreted as yielding the ground- 
state configuration that corresponds to the jj-coupled Slater determinant formed of six j = 
5/2 orbitals. Further, taking into account the finite (« 2 eV) width of the /-band and the 
occupation U/ = 5.44, we can interpret the calculated ground state as a partially localized / 6 
manifold hybridized with a broad valence band (see Fig. 3(a)). 

An important criterion for assessment of results of electronic structure calculations is given 
by electron spectroscopies. Although some techniques, as BIS, have not been applied on Pu 
yet, there exists high-resolution Ultraviolet Photoelectron Spectroscopy (UPS) data, mapping 
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the electronic structure down to about 10 eV below Ep [21,23-25], while X-ray photoelectron 
spectroscopy [21,24,25] gives information on the screening of a deep core-hole by conduction 
electrons. 

Valence-band spectra of S— Pu and handful of other Pu-based systems studied so far exhibit 
invariably three narrow features, one at Ep or in its close vicinity, the other two at 0.5 and 0.85 
eV below Ep, which are called A, B, and C, respectively. Some systems exhibit also another, 
broader, separate feature D, the energy of which varies between compounds, and which is 
not noticeable for S— Pu. It was argued [24] that the intensity of D reflects the degree of 5/ 
localization. For example, the peaks A-C are absent in the 5/ localized compound PuSb [26], 
and all the 5/ spectral intensity is concentrated in the broad peak D. 

Valence-band photoelectron spectra reflect to some extent the density of occupied states 
in the ground state (DOS). A comparison with the bulk experimental UPS [21,24] (see Fig. 
3(c)) shows that LDA+U places the 5/ manifold approximately 0.9 eV below the Fermi energy 
(Ep) in accordance with the experimental C peak position, while it does not resolve correctly 
the experimental features A and B at the Ep edge. 

As discussed in Ref. [22], the sharp features A-C are not specific for individual compounds, 
but they appear at invariable energies over a class of diverse Pu systems. It is therefore 
implausible that they could be related to any ground-state one electron states. Dismissing 
the peaks A-C from the comparison, as we assume them either to be due to the final-state 
effects, which could be reconciled with the 5/ 5 final-state multiplet [27], or due to a general 
many body phenomenon (analogy of a Kondo effect) [28]), we are left at a difficult situation, 
as no other well-defined feature is seen in the UPS spectra. But for ultrathin layers, with 
thickness varying down to about one monolayer (see Fig. 3(c)), a broad, but well defined 
peak D, centered at 1.6 eV binding energy, concentrates most of the 5/ spectral intensity, 
whereas the peaks A-C (still at their energies) become barely noticeable. 

Therefore we undertook analogous AMF-LDA+U calculations for free-standing Pu slabs 
(where the number of Pu-layers varies from 1 to 5) with the lattice parameter identical to 
that of (5— Pu. For the one monolayer (1ML) (see Fig. 3(d)), the 5/ states form a very narrow 
band at 1.5 eV, which is in an excellent agreement with the position of the D-peak (cf. Fig 
3(c)). An intermediate situation, Pu trilayer (3ML), demonstrates that increasing thickness 
leads to the band broadening, but also to a shift towards Ep. With a further increase of the 
Pu slab thickness up to five monolayers (5ML), the band broadening continues and the DOS 
becomes very similar to the bulk (cf. Fig. 3(a)). This broader 5/ band is in fact underlying 
the peak C, which explains why it is not observed. Moreover, the fact that in all four cases 
(bulk, 1, 3, and 5 monolayers) the non-magnetic state (L = , S = 0) persists, confirms the 
robustness of this solution. 

It is also to mention that both LSDA and FLL-LDA+U magnetic solutions yield the DOS 
which does not correspond to the experimental UPS spectra shown in Fig. 3(c). In the LSDA 
case, the 5/-band has a very sharp peak structure in the region from —1.2 eV to 0.5 eV 
around Ep. For the FLL-LDA+U, the 5/- band becomes localized with the 5/ manifold well 
separated below the bottom of the valence band in the energy region from 2.5 to 4.5 eV below 
the Ep, again in a contradiction to experimental data. 

To summarize, we have shown that the around-the-mean-field version of the LDA+U 
method yields non-magnetic ground state of 5— Pu and, at the same time, it gives a good 
agreement with experimental data, namely, with the photoemission spectra, equilibrium vol- 
ume, and bulk modulus. The fact that the non-magnetic character is not due to a mutual 
cancellation of spin and orbital parts of the moment, but is due to S — and L — , is par- 
ticularly corroborated by very recent neutron scattering experiments [5]. We are aware that 
our approach does not account for dynamical effects, which may become very important, as 
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suggested in Ref. [9]. Also, as we focused here on the magnetic, structural and spectroscopic 
properties of the 5— Pu only, we do not discuss the a— to Pu phase transition, leaving it 
for further consideration. 
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